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Abstract 

A constitutive model is derived for the time-dependent response of particle-reinforced elastomers 
at finite strains. An amorphous rubbery polymer is treated as a network of long chains linked by 
permanent junctions (chemical crosslinks, entanglements and filler particles). A strand between two 
neighboring junctions is thought of as a sequence of mers whose motion is restricted to some tube by 
surrounding macromolecules. Unlike the conventional approach that presumes the cross-section of 
the tube to be constant, we postulate that its radius strongly depends on the longitudinal coordinate. 
This implies that a strand may be modeled as a sequence of segments whose thermal motion is 
totally frozen by the environment (bottle-neck points of the tube) bridged by threads of mers which 
go through all possible configurations during the characteristic time of a test. Thermal fluctuations 
affect the tube's radius, which results in freezing and activation of regions of high molecular mobility 
(RHMs) . The viscoelastic response of an elastomer is associated with thermally activated changes in 
the number of RHMs in strands. Stress-strain relations for a rubbery polymer at finite strains and 
kinetic equations for the concentrations of RHMs are developed by using the laws of thermodynamics. 
At small strains these relations are reduced to the conventional integral constitutive equation in linear 
viscoelasticity with a novel scaling law for relaxation times. The governing equation is determined 
by 5 adjustable parameters which are found by fitting experimental data in tensile dynamic tests on 
a carbon black filled natural rubber vulcanizate. 

Key-words: Particle-reinforced elastomers, Viscoelasticity, Tube concept, Thermal fluctuations, Dy- 
namic tests 

1 Introduction 

This study is concerned with the viscoelastic behavior of unfilled and particle-reinforced rubbers in 
isothermal tests. The time-dependent response of elastomers and rubbery polymers has attracted essential 
attantion in the past decade, which may be explained by numerous applications of rubber-like materials 
in industry (vehicle tires, shock absorbers, earthquake bearings, seals, flexible joints, solid propellants, 
etc.). Despite noticeable progress in this field (see, e.g., Govenjee and Simo, 1992; Ozupek and Becker, 
1992; Gonzonas, 1993; Johnson and Stacer, 1993; Le Tallec et al., 1993; Zdunek, 1993; Hausler and Sayir, 
1995; Johnson et al., 1995; Aksel and Hiibner, 1996; Holzapfel and Simo, 1996; Lion, 1996, 1997, 1998; 
Spathis, 1997; Bergstrom and Boyce, 1998; Ha and Schapery, 1998; Reese and Govinjee, 1998; Septanika 
and Ernst, 1998a,b; Ernst and Septanika, 1999; Mieche and Keck, 2000; Wu and Liechti, 2000; Haupt 
and Sedlan, 2001), some important questions remain, however, obscure. 

Numerous observations in uniaxial and biaxial tests evidence that commercial rubbery vulcanizates 
reinforced with carbon black and silica, as well as solid propellants filled with micrometer-size hard par- 
ticles demonstrate a pronounced time-dependent response which may be associate with their viscoelastic 
properties. For examples, we refer to recent studies by Hausler and Sayir (1995), Lion (1996-1998), 
Bergstrom and Boyce (1998), Septanika and Ernst (1998b), Mieche and Keck (2000), Wu and Liechti 
(2000), and Haupt and Sedlan (2001). An unresolved problem remains about the mechanism for stress 
relaxation at the molecular lovel. 
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To answer this question, Aksel and Hiibner (1996) carried our tensile relaxation tests on crosslinked 
polybutadiene rubber reinforced with glass beads at ambient temperature with various volume fractions 
of filler (from to 45 %). Obvervations evidenced that there was no relaxation in the unfilled rubber, 
and the level of relaxed stress increased with the content of reinforcement. They suggested that slippage 
of polymeric chains along the filler surface may serve as a mechanism for stress relaxation in the linear 
regime, whereas detachment of chains from the surfaces of beads and creation of valuoles characterized 
the nonlinear viscoelasticity of rubber. 

Liang et al. (1999) reported experimental data in dynamic tensile tests (with frequency 1 Hz) on 
ternary composite of polypropylene and EPDM (ethylene propylene diene monomers) elastomer filled 
with glass beads. They demonstrated that (i) treatment of beads by silane coupling agent did not affect 
the mechanical damping of specimens (measured in terms of the loss modulus) and (ii) the loss modulus 
of the compound with untreated beads decreased with the concentration of filler [in contrast with the 
observations by Aksel and Hiibner (1996)]. 

These conclusions were recently questioned by Leblanc and Cartault (2001), who measured dynamical 
shear response of uncured styrene-butadiene rubber reinforced with various amounts of carbon black 
and silica and found that the compounds exhibited practically no viscoelasticity or rather limited time- 
dependent response. This implies that crosslinks play the key role in viscoelasticity of rubber, while the 
effect of filler particles is relatively weak. 

Clarke et al. (2000) carried out tensile relaxation tests on two unfilled elastomers (natural polyisoprene 
rubber with polysulphur crosslinks and synthetic polyacrylate rubber with carbon-carbon and carbon 
oxygen crosslinks) at the temperature T — 80 °C. At a relatively small strain (e = 0.05) the natural 
rubber demonstrated significant mechanical relaxation (which is observed during at least 10 6 s), whereas 
the synthetic elastomer revealed no relaxation. Clarke et al. (2000) associated slow relaxation of natural 
rubber with isomerisation of polysulphur crosslinks, but could not explain stress relaxation in natural 
rubber observed in short-term mechanical tests, as well as its absence in polyacrylate elastomer. 

Two approaches may be mentioned to the description of the time-dependent response in rubbery 
polymers. The first was proposed about half a century ago by Green and Tobolsky (1946) and, afterwards, 
revised by Yamamoto (1956), Lodge (1968) and Tanaka and Edwards (1992), to mention a few, and 
applied to describe the viscoelastic behavior of rubbery polymers by Drozdov (1997) and Ernst and 
Septanika (1999). According to the concept of transient networks, a rubbery polymer is modeled as a 
network of long chains bridged by temporary junctions (which are associated with chemical and physical 
crosslinks and entanglements). A strand whose ends are linked to neighboring junctions is treated as an 
active one. Snapping of an end of a strand from a junction (breakage of a crosslink or disentanglement 
of a chain) is thought of as its breakage (transition from the active state of the strand to its dangling 
state). When a dangling chain captures a junction (which reflects creation of a new physical crosslink or 
entanglement), a new active strand merges with the network. Breakage and reformation of active strands 
arc treated as thermally activated processes: attachment and detachment events occur at random times 
as they are agitated by thermal fluctuations. A disadvantage of this model is that it does not explain 
observations by Leblanc and Cartault (2001): the theory of transient networks implies an increase in 
the strength of relaxation in reinforced rubbery compounds independently of the degree of crosslinking, 
because particles of filler play the role of extra junctions in the network. It is worth mentioning two 
shortcomings of this theory from the theoretical standpoint: (i) it associates viscoelasticity with the 
micro-motion of strands, which contradicts the conventional opinion that relaxation of stresses reflects 
local segmental motion in the chain backbone (Sperling, 1996), and (ii) this approach is grounded on the 
hypothesis that the stress-free state of a reformed strand coincides with the deformed state of the network 
at the instant of their merging. The latter is tantamount to implicit ascribing some elastic features to 
long chains [see a recent discussion of this issue for microemulsions with telechelic polymers in Filali et al. 
(2001)], whereas the conventional theory of rubber elasticity disregards the contribution of mechanical 
energy into the free energy of chains compared to that of configurational entropy (Treloar, 1975). 

The other concept was recently suggested by Clarke et al. (2000) who associated the viscoelas- 
tic behavior of elastomers with mechanically induced rearrangement (isomerisation) of crosslinks. The 
shortcomings of this approach are that (i) it contradicts experimental data by Aksel and Hiibner (1996) 
(because it disregards the effect of filler) and (ii) this model does not explain the presence of differ- 
ent relaxation modes observed in a number of mechanical tests (because crosslinks are identical, their 
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rearrangement results in a unique relaxation time inversely proportional to the characteristic time of 
rearrangement), and (iii) it fails to describe the viscoelastic behavior of uncrosslinked polymers (uncured 
rubbers or thermoplastics above the glass transition temperature). 

The objective of this study is three-fold: (i) to develop a micromechanical model which does not 
contradict the above experimental data, (ii) to find adjustable parameters in constitutive equations by 
fitting observations in dynamic tests with a wide range of frequencies (3 decades) at various temperatures 
(from —20 to 100 °C), and (iii) to demonstrate that the material constants are altered with temperature 
in a physically plausible way. By no means this model is treated as an alternative for existing theories, 
but rather as an accompanying concept. This is a reason why we focus on observations in short-term 
tests (the frequency, /, changes from 10 _1 to 10 2 Hz) with relatively small strains (less than 10 %), 
by assuming that (i) the theory of transient networks with mechanically induced rates for breakage and 
reformation of active strands can correctly predict experimental data in the region of larger strains, where 
disentanglement of chains and dewetting of the host polymer from the filler surfaces become essential, and 
(ii) the theory of junctions' rearrangement can be adequately applied to the analysis of the time-dependent 
response in long-term tests (with the characteristic time of the order of 10 6 s). 

At first glance, approximation of experimental data in isothermal dynamic tests with small strains is 
an oversimplified task. Assuming the stress, a(t), at an arbitrary time t > to depend linearly on the 
strain history, e(r), r G [0, t], and referring to the Risz theorem for the presentation of a linear functional, 
one can write 



Jo J 

where E is an elastic modulus and R(t) is a relaxation kernel. Expanding the function R(t) into a 
truncated Prony series, 



one finds that given a number of terms, N\, and a series of relaxation times, tn, the parameters E and £jv 
may be uniquely determined by matching the storage and loss moduli with the help of the least-squares 
algorithm. This approach perfectly works when we aim to approximate observations at some temperature 
and amplitude of oscillations. In a more complicated case which is addressed here, when observations are 
available at several temperatures and amplitudes of strains (the response of particle-reinforced rubber 
vulcanizates is extremely sensitive to the amplitude of oscillations even for strains of the order of several 
per cent) and the objective is to evaluate the effect of thcrmomechanical factors on material parameters, 
phcnomcnological equations (1) and (2) totally fail, because they provide no information about changes 
in the relaxation times, tjv, with temperature and strain. 

To shed some light on the effects of temperature and strain on the relaxation spectrum, we develop 
a model for the time-dependent response of elastomers at finite strains (i) that is reduced to Eq. (1) in 
the case of small strains, and (ii) which expresses the quantities £jv and tn in Eq. (2) in terms of some 
micromechanical parameters, whose dependence on thermomechanical factors appears to be physically 
plausible. This model contains only five material constants which are determined by fitting observations. 

According to the statistical theory of macromolecules, a rubbery polymer is thought of as a network 
of long chains linked by permanent junctions (crosslinks, entanglements and van der Waals forces). A 
strand between two neighboring junctions is thought of as a set of large number of mers connected in 
sequel. The classical theory of rubber elasticity treats an individual strand as an object independent of its 
surrounding, which implies that the configurational entropy (associated with the number of configurations 
available for statistically independent segments in the strand) provides the main contribution into the free 
energy per strand. It was recognized four decades ago that micro-motion of segments is confined to the 
neighborhood of some "mean configuration" by surrounding macromolecules which suppress their thermal 
fluctuations. This confinement is traditionally accounted for in the framework of the tube concept which 
presumes that a chain can freely move only within some tube created by neighboring macromolecules 
(Doi and Edwards, 1986). For a detailed discussion of constrained fluctuations, the reader is referred to 
Everaers (1998). 

The tube is conventionally treated as a curvilinear cylinder with a constant cross-sectional area. The 
main hypothesis of the present study is that small-scale local heterogeneities dramatically affect micro- 
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motion of a strand, which implies that the radius of the cylinder strongly depends on the longitudinal 
coordinate. This means that the surface of the cylinder is extremely wavy, and the tube may be treated as 
a sequence of regions with high mobility of segments bridged by bottle-neck regions, where the molecular 
mobility is strongly suppressed. The sketch depicted in Figure 1 demonstrates the position of the tube and 
the location of regions of high mobility (RHMs) linked by regions of suppressed mobility at the current 
time t. Thermal fluctuations of surrounding macromolccules change the positions of RHMs noticeably 
even after a small interval of time, At, which implies that the current properties of a strand may be 
described in terms of the statistical distribution of RHMs only. 

To make the model tractable, we introduce the distribution, pj^, of strands with various numbers of 
RHMs, N, in an equilibrated sample (in the stress- free state at time t = 0). Application of external 
loads activates the network, which implies that the concentrations of RHMs in strands are altered. At an 
arbitrary instant t > 0, the state of a strand (which initially has N RHMs) is described by the current 
number of regions of high mobility, N a (t). Suppression of RHMs by surrounding molecules, as well as 
activation of frozen segments (agitated by thermal fluctuations occuring at random times) are assumed 
to serve as a basic mechanism for the viscoelastic response of an elastomer at the micro-level. 

The study is organized as follows. Deformation of a strand is described in Section 2. The mechanical 
energy of a network is determined in Section 3. Stress-strain relations at finite strains and a kinetic 
equation for the function N a (t) are derived in Section 4 by using the laws of thermodynamics. The 
constitutive equations are simplified for the case of small strains in Section 5. Uniaxial tension of a 
specimen is analyzed in Section 6. Results of numerical simulation are compared with observations in 
dynamic tests on a carbon black filled natural rubber vulcanizate in Section 7. The experiments exposed 
in the introductory section are shortly discussed in terms of the model in Section 8. Some concluding 
remarks arc formulated in Section 9. 



2 Deformation of a strand 

We consider a strand linking two junctions, A\ and A2, which contains N regions of high molecular 
mobility in the stress-free state at the initial instant t = and N a (t) RHMs at time t > 0. Introducing 
the dimensionless increment of the number of RHMs, 

»«) = ^« 

we find that 

N a (t) = N[l - n(t)}. (3) 

It is worth noting that the function n(t) may accept both positive (N a (t) < N) and negative (N a (t) < N) 
values. 

Let e be a strain from the stress- free state of a strand to its deformed state. Assuming that suppressed 
regions arc not deformed, we find that the strain e equals the sum of strains for active RHMs, 

e = N a e a . 

This equality together with Eq. (3) implies that 

ea ^ NlT^-y (4) 

The mechanical energy of suppressed regions (totally frozen segments of a strand) is assumed to vanish. 
Active RHMs are modeled as linear elastic solids with the strain energy 

T^a, (5) 

where the rigidity per bond [i has the dimension of energy. Two remarks are necessary regarding Eq. 
(5). First, this relation presumes small deformations of RHMs. This hypothesis was confirmed by recent 
molecular dynamic simulations (Bcrgstrom and Boyce, 2001) which demonstrated that axial elongation 
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of a specimen with macro-strains up to 100 % implied micro-strains per node which did not exceed 5 
% (for chains containing 200 atoms, which is in quantitative agreement with data obtained by fitting 
observations in our study, see Section 7). Second, at small strains, the contributions of mechanical and 
statistical elasticity into the free energy of a polymer cannot be distinguished (Ferry, 1980), which implies 
that n = n(T) is the rigidity associated with both (entropic or mechanical) components of the potential 
energy of deformation. 

The strain energy of a strand, wn, is determined as the sum of the strain energies for individual 
RHMs. It follows from Eqs. (4) and (5) that 

"»V = l^l = W^y (•> 

To express the strain in a strand, e, in terms of the Cauchy deformation tensor for the network, C, we 
suppose that in the stress-free state the strand has a small end-to-end length S n and it is directed along 
the unit guiding vector I, 

f°-fi= 6 N l, 

where r° and are initial radius vectors of the end points, A\ and A-2,. At some instant t > 0, the 
junctions occupy points with the radius vectors 

n (t) = fl + u(t, f? ) , r 2 (t) = f 2 ° + u(t, f° ) , 

where u(t, f) is the displacement vector at point f for transition from the stress-free state to the deformed 
state at time t. The current position of the strand is determined by the end-to-end vector 

R(t) = f 2 (t) - fi(t) = 6 N l + [u(t, ff + d N l) - u(t, fl)}. 

Neglecting terms beyond the first order of smallness compared to Sjy, we find that 

R(t) = S N l- [J + V°ti(t,f?)], 

where V° is the gradient operator in the stress-free state, I is the unit tensor and the dot stands for inner 
product. In terms of the radius vector in the actual state, f, this equality reads 

R(t) = S N L V°f(t) = 5 N [V°f{t)] T ■ I, (7) 

where T stands for transpose. The end-to-end length of the strand, ds, is given by 

ds 2 = RR. 

It follows from this equality and Eq. (7) that 

ds 2 (t,l) = S 2 N l-C(t)-l, (8) 

where 

CW=V°f(i)-[V r-(*)] T (9 ) 

is the Cauchy deformation tensor for transition from the stress-free state of the network to its deformed 
state at time t. 

The extension ratio for a strand, A, equals the ratio of the current end-to-end length, ds, to that 
of the strand in an "activated" stress- free state, ds° . The latter state is a state of the strand which is 
suddenly unloaded, but for which the number of RHMs coincides with its current value, iV a (i). It differs 
from the equilibrium stress- free state of the strand, where the number of RHMs equals N. The difference 
between the end-to-end lengths of the strand in the activated state, ds°(t), and in the equilibrium state, 
ds(0) = Sn, determines the end-to-end elongation driven by changes in the number of RHMs. We assume 
that the transformation-induced end-to-end elongation of a chain is proportional to the number of frozen 
RHMs, N - N a (i), which implies that 

ds°(t,l) = S N + Nn(t,T)6° , (10) 
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where 5° is an increment of the end-to-end length driven by transition of an individual RHM. It follows 
from Eqs. (8) and (10) that 

ds(t,l) _ [l-C(t)-l]i 



A(t,i) 



where 



ds°(t,l) l + t] N n(t,l)' 
NS° 



(11) 



VN 



Adopting the conventional hypothesis in rubber elasticity (which is tantamount to the assumption about 
statistical independence of RHMs, see Treloar, 1975) that 

s N = sVn, 

where the coefficient 5 is independent of N, we arrive at the formula 

. 5° 



V s ■ 



■q N = fjvN, 

The strain in a strand is determined by the Hencky formula 

e = In A. 

Substitution of expression (11) into this relation implies that 



(12) 



e(t,0 = ^ln l-C(t)-l 



In 



l + r] N n(t,l) 



(13) 



Equation (13) establishes correspondence between the micro-deformation of a strand which is described 
by the strain, e, and the macro-deformation of a network which is determined by the Cauchy deformation 
tensor, C. The parameter rjjy in Eq. (13) may serve as a measure of coiling for strands: for strongly 
coiled chains, transformation- induced elongation of a strand is negligible compared to its length and 77 jv 
is close to zero, whereas for uncoiled chains t]m has the order of unity. 



3 Strain energy density of a network 

In the stress-free state a rubbery polymer is modeled as a network of strands with various numbers 
of RHMs (N — 1, 2, . . . , N\), where N\ is the maximal number of regions of high mobility per strand. 
Denote by En the concentration (per unit mass) of strands with TV RHMs (in thermal equilibrium in the 
stress-free state) and by 

N=l 

the total number of strands per unit mass. The distribution of strands is characterized by the function 

Pn = — ■ 

Referring to the random energy model (see, e.g., Dyre, 1995), we accept the quasi-Gaussian ansatz for 
the probability density pn, 

(N-N) 2 - 



Pn = Aexp 



2Y, 2 



(14) 



where N and S are positive constants, and the parameter A is determined by the condition 



£f>iv = 1. 

N=l 
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When N » S, we can identify N with the average number of RHMs in a strand and £ with the standard 
deviation of the number of RHMs. 

We accept the conventional assumption that the excluded-volume effect and other multi-chain effects 
are screened for an individual chain by surrounding macromolecules (Everaers, 1998). This implies that 
the energy of interaction between strands may be neglected (under the hypothesis of incompressibility) 
and the strain energy of the network may be determined as the sum of the potential energies of deformation 
for individual strands. Assuming the distribution of strands with various guiding vectors to be isotropic, 
we find the concentration (per unit mass) of strands with N RHMs and guiding vector I, 

Xn(1) = -j—EpN siwddfldip, (15) 
Air 

where and ip are Euler's angles that determine the position of the unit vector [. Multiplying the number 
of strands, X^(l), by their mechanical energy, wn, summing the results for various ./V and I, and using 
Eqs. (6) and (15), we find the strain energy of the network 

8 ^ hi N Jo Jo 1-«M,^) 1 ' 

Our purpose now is to find the derivative of the function W with respect to time. Because the parameters 
/i, S and pn are, in general, functions of temperature, T, the derivatives of these quantities should be 
accounted for in the calculation of dW/dt for an arbitrary thermo-mechanical process. To simplify the 
analysis, we confine ourselves to small changes in temperature and presume that the concentrations of 
strands, Ejv, with various numbers of RHMs, TV, are affected by temperature rather weakly. This implies 
that only the derivative of the average rigidity per bond, /i, is taken into account in the formula for the 
derivative of W. Differentiation of Eq. (16) results in 



- ek r a r e M^) 

dt ^ki N Jo Jo l-n{t,#,<pY 



' U -" T "' ■^(t,0,<p) S m0d0 + J o (t)^(t) + .hil). (17) 



where 



" jfc N Jo Jo l-n{t,0,<p) 



N= 

/iS A p N f 271 , f n e 2 (t,d,(p) dn, „ . . „,„ 



(18) 



The function Jo(t) in Eq. (18) determines the increment of the strain energy of a network associated 
with the dependence of the average rigidity per RHM, /i, on temperature, whereas the function Ji(t) 
characterizes the increment of the potential energy of deformation driven by suppresion and thermal 
activation of RHMs. Differentiation of Eq. (13) results in 



de, 1 



dC 



-tt,i) = ^^ — =■ i.^L( t )-i] —, T ^(t,n. (19) 

dV'' 2[l-C(t)-l}\- dt K> J l + m n(t,l) dt K ' v ; 



It follows from Eq. (9) that 

dC 



where 



<u (t) = V°«(t) • [V°f(i)] T + V°f (*) • [V°«(t)] T , 



m = Tt (t) 



is the velocity vector for the network. Bearing in mind that (Drozdov, 1996) 

V°u(t) = V°f(t) • V(t)v(t), 
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where Vft) is the gradient operator in the deformed state at time t, we obtain 

dC , 



dt 



where 



ft) = 2V°rft) • D(t) ■ [V°fft)] 1 , 



l 



D(t) = -[W(t)v(t) + (W(t)v(t)) 



is the rate-of-strain tensor for the network. Equation (20) implies that 



(20) 



where 



dC 

I . — ft) -1 = 21- V°fft) • D(t) • [V°f ft)] T • Z = 2Fft, : £>(t), 



Fft, i) = [V°f ft)] T • (/ ® I) • V°f ft) 



(21) 



(22) 



is the generalized Finger tensor, the colon stands for convolution and ® denotes tensor product. Substi- 
tution of Eqs. (19) and (21) into Eq. (17) results in the formula 



^ft) = Tft):Dft) + J ft)^ft)-J 2 ft), 



where 



2tt ,-t: 



dip 



4tt ^ N J Jo 1 - nft, i?, f . Cft) . [ 



sin #cfa?, 



M ~ ^ PJV , r mm,- 



2-TT />7T 



h N (t 7 d 7 ip) dn 



h N (t,i),(p) = eft,i?,p) 



27?^ [1- raft, t?, y)] 
1 + Jjjvnft, i?, V 5 ) 



eft 



(23) 



(24) 



4 Constitutive equations 

Observations evidence that the growth of temperature driven by periodic straining with small amplitudes 
(less than 10 %) is quite modest [it does not exceed 10 to 20 K for frequencies, /, in the range from 10 _1 to 
10 2 Hz, see Kar and Bhowmick (1997) and Lion (1998)]. This implies that the absolute temperature, T, 
remains rather close to its reference value, To, and thermal expansion of the network may be disregarded. 
Under the hypothesis about affine deformation (which states that the rate-of-strain tensor for deformation 
of a specimen coincides with the rate-of-strain tensor for the network, D) , the first law of thermodynamics 
for an incompressible network reads (Coleman and Gurtin, 1997) 

- = -(* d :D-V.,) + r, (25) 

where p is mass density, is the deviatoric component of the Cauchy stress tensor a, q is the heat flux 
vector, $ is the internal energy and r is the heat supply per unit mass. At relatively small increments 
of temperature, the density of a polymer, p, is treated as a constant. The Clausius-Duhem inequality is 
given by (Coleman and Gurtin, 1967) 

d_Q = dS ,|n _pr > 

H dt 1 dt \TJ T ~ ' K ' 

where S is the entropy and Q is the entropy production per unit mass. Bearing in mind that 

V-(|) = lv- 9 --^--VT 
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and excluding the term V • q from Eqs. (25) and (26), we find that 

The internal energy per unit mass, is connected with the free (Helmholtz) energy per unit mass, 
by the equality 

$ = f + ST. (28) 

It follows from Eqs. (27) and (28) that 

dQ dT d* 1/ . 1.^. „ 

T ^=- S ^-^ + - p { ad:D -T q - VT )^°- (29) 

It is assumed that the free energy, 'J, consists of the following four terms: 

1. the free energy, ^oq, of the network in the equilibrium stress- free state at the reference temperature 
To (when the numbers of RHMs in a strand coincides with its initial value, TV); 

2. the free energy of thermal motion, ^ t h, which is determined by the conventional formula 

*th = (c-So)(T-To)-cTln^, 

where c is the specific heat (at relatively small increments of temperature, c is treated as a constant) 
and So is the entropy per unit mass in the stress- free state at the reference temperature To; 

3. the potential energy of deformation for the network, W, given by Eq. (16); 

4. the increment of free energy, A 'J, induced by suppression and thermal activation of RHMs. 

The latter term is associated with the excess free energy of a strand with the current number of RHMs, 
iV a , over that for the strand with the equilibrium number of RHMs, N. As a first approximation, the 
increment of the free energy for a strand with N RHMs, A^jv, is given by 

AV>Ar = i/?/i(^a-A0 2 , (30) 

where j3 is a dimcnsionless parameter (which, in general, depends on temperature), and the average 
rigidity per RHM, ^, is introduced to ensure the correct dimension of the function Aipjy. Summing the 
increments of free energy, Eq. (30), for strands with various numbers of RHMs, N, and various guiding 
vectors, I, and using Eqs. (3) and (15), we find that 

BuF Nl f 2lT r 

= ^5= N 2 p N / dip n 2 (i,i?,<^)sintf*?. 
8^ ^ Jo Jo 

The free energy per unit mass of the network, reads 

T Bu" Nl f 27r 

* = * cq +(c-S*o)(T-T )- C Tln— + W+^ ^ N 2 p N / dtp \ n 2 sinzW. (31) 

It follows from Eqs. (23), (24) and (31) that the derivative of the free energy per unit mass, with 
respect to time, t, is given by 

§-t:A + (^-*-ch£)§-J, (32) 
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where 

Nl r 27T 



AS(t) = Jo (t) + f - V JV=W ["dtp [" n 2 sin tfdtf, 

8tt ^ „/ J 

87r ^ N Jo Jo [l-n{t,d,<p)\ 2 dt 



H N {t,0,<p) = e(t, 0, cp) \ 2llN ^ H ^ p ^3 - eft, fl, y) } - 2/37V 3 n(t, fl, y) [1 - n(t, fl, y)] 2 . (33) 

L 1 + VNn(t, v, ip) J 

Combining Eqs. (29) and (32), we find that 

T j4(»'-' t ) :i -( s -*- t £ + as ) § - w» -* t + ^ »■ < 34 » 

Applying the conventional reasoning to Eq. (34), see Coleman and Gurtin (1967), we find that the 
expressions in braces vanish. This implies the formula for the entropy per unit mass, 

S = S + cln^-AS (35) 
Jo 

and the relationship for the Cauchy stress tensor 

where P(t) is pressure. It follows from Eq. (35) that the quantity AS* [given by Eqs. (18) and (33)] char- 
acterizes changes in the entropy per unit mass, S, driven by deformation of the network and suppression- 
activation of RHMs. Substitution of Eq. (22) into Eq. (36) results in the constitutive equation 

&(t) = -P(t)I + [V°f(t)] T • B(t) • V°f(t) (37) 

with 

PH~ ^ p N f 27T J r e(t,d,(p) I® I 



B(t) = ^y^L d(p ^ _ sin ^. (38 ) 

4n ^i N Jo Jo l-n(t,#,<p)l.C(t)-l V ' 

Combining Eqs. (34) to (36), we obtain 

The function J(t) determines the rate of entropy production (per unit mass) induced by changes in 
the numbers of RHMs, whereas the last term on the right-hand side of Eq. (39) characterizes entropy 
production (per unit mass and unit time) driven by thermal conductivity. The function Q is assumed to 
be presented in the form 

Q(t) = Qth(t) + AQ(t), (40) 

where Q t h describes the entropy production induced by thermal diffusivity and AQ determines the entropy 
production associated with the suppression-activation mechanism. We accept the conventional formula 
for the derivative of the function Qth, 

dt () ~ P T* ' 1 > 

where n > is a constant thermal diffusivity, and (VT) 2 = VT • VT. It follows from Eqs. (39) to (41) 
that 



T dAQ = J 



dt pT 



-L(> + kVt) -VT>0. (42) 
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Bearing in mind that this inequality holds for an arbitrary thermo-mechanical process, we find that the 
expression in braces vanishes, which results in the Fourier law for the heat flux vector 

q{t) = -nV{t)T(t). (43) 

Equations (42) and (43) yield 

T ^ir = J ^ °- (44) 

We postulate that the function AQ is given by 

AQ(t) = — Y~ N dcp AQ N (t,0,<p)fan#d#, (45) 
4tt ^ Jo J 

where AQtv(£, I) characterizes the energy dissipation per strand containing N RHMs and directed along 
vector I. To construct the functions AQn, we refer to the conventional expression for the entropy 
production Q in an incompressible viscous fluid 

T^§ = (D: D, (46) 

where £ stands for viscosity. Equation (46) means that the rate of entropy production per unit mass is 
proportional to the square of the relative velocity (which is characterized by the rate-of-strain tensor D). 
By analogy with this relation, we postulate that the expression 

T(t)^M,^) 

is proportional to the square of the relative rate of changes in the number of RHMs. According to Eq. 
(3), this rate is given by 

_ _l_cWa = _J_dn 

N a dt 1-ndt' [ ' 

With reference to Clarke et al. (2000), we assume that energy can dissipate at junction points only, which 
implies that Eq. (47) should be multiplied by the probability that the activated (suppressed) segment 
is located at the end of a strand (this probability equals 2/./V for a strand containing N RHMs). As a 
result, we obtain 

T(t)^L(t, tf.p) = ^ [l - n(t, 0, <p)] 2 \*£(t, 0, <pj\ 2 , (48) 



dt 



where k is a positive coefficient independent of the number of RHMs in a strand. Substituting expressions 
(45) and (48) into Eq. (44) and equating coefficients at the derivative dn/dt on the left and right hand 
sides on the resulting equality, we arrive at the kinetic equation for the function n(t, 0, <p), 

It follows from this formula and Eq. (33) that 

^(t, 0, if) = a 1N e(t, 0, l){ 2 ^~ n y^ )] - e(t, 0,?)} - ^n(t, 0, V )[l - n(t, 0, ^)] 2 , (49) 



where 



r4 M M/3 



aiN = aiN, a 2 N = a 2 N , ai = — , a 2 = — . (50) 
It follows from Eq. (3) that the initial condition for Eq. (49) reads 

n{t,d,<p) =0. (51) 
t=o 

Relations (13), (37), (38), (49) and (51) provide a set of constitutive equations for a rubbery polymer at 
finite strains. The behavior of the coefficients (50) in kinetic equation (49) is unusual. The parameter aiN 
linearly increases with N , which means that suppression of RHMs occurs more rapidly in long strands 
than in short ones. On the other hand, the quantity a 2 N grows as the fourth power of N, which implies 
that the concentration of RHMs in long strands remains sufficiently large. 
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5 Constitutive equations at small strains 



We now simplify the constitutive equations for the case of small strains. According to Eq. (49), the 
condition \e\ <C 1 ensures small changes in the concentrations of RHMs, n«l. Introducing the strain 
tensor e by the conventional formula 



we find from Eq. (13) that 



e(t,l) = -ln|l + 2Z-e(t) -l 



- In 



l + r) N n(t,l) 



Neglecting terms beyond the first order of smallness compared to |e| and n, we obtain 

e(t, I) — I ■ e(t) ■ I — t]Nn(t 7 1). 



(52) 



It follows from Eqs. (49) and (52) that with the required level of accuracy n satisfies the differential 
equation 



^(i, I) = 2ai N r) N l ■ e(t) ■ I - (2ai N r] 2 N + a 2N )n{t, I). 
The solution of Eq. (53) with initial condition (51) reads 

n(t,l) = 2aiNi]N / exp -(2aiNVN + a 2N)(t - r) I ■ e(r) • I dr. 
Jo 1 J 

This equality together with Eq. (52) implies that 

e(t, I) = I - \ e(t) - 2a lN Tf N / exp -~(2ai N ri 2 N + a 2 jv)(i - r) e(r) dr \ ■ I. 



(53) 



(54) 



Substituting Eq. (52) into Eqs. (37) and (38) and neglecting terms beyond the first order of smallness, 
we find that 



JV=1 



dtp / {l®l)sm i ddd 



L 



I ■ \ e(t) - 2a 1N r] 2 N I exp -(2a 1N rj 2 N + a 2 N){t - r) e(r) dr\ ■ I. 



(55) 



To calculate the integral, we introduce a Cartesian coordinate frame {xi} with base vectors e, which are 
directed along the eigenvectors of the symmetric tensor e and obtain 



e(t) = e 1 (t)e 1 e 1 + e 2 (t)e 2 e 2 + e 3 (t)e 3 e 3 , 
where is the ith eigenvalue of e. The unit vector I is given by 

I = cos i?ei + sin $ (cos ipe 2 + sin ipe 3 ) , 
It follows from Eqs. (56) and (57) that 

I ■ e(t) ■ I — ei (t) cos 2 $ + e 2 (t) cos 2 ip + e 3 (t) sin 2 ip 

Equation (57) implies that the tensor F® I is determined by the matrix 



(56) 



(57) 



sin 



I (g) I 



cos 



sin 1? cos i? cos <p> sin ■& cos § sin (p 



sin $ cos •& cos ip sin i9 cos tp sin d sin ip cos ip 
sin i? cos sin y> sin 2 $ sin <p cos sin 2 $ sin 2 <p 
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After simple algebra we find that 



271 



4^1 diP 



(l®l)sm$dtf 1-1-1 = — 
15 



3ei + e 2 + £3 

ei + 3e 2 + £3 














£i + £2 + 3f 3 



This equality together with the incompressibility condition, 

£l + £2 + £3 = 0, 

results in the formula 



2tt 



to 1 ^ 



(Z (8) sin Mi? 1-1-1= —e. 
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Substituting this expression into Eq. (55), we arrive at the stress-strain relation 

s{t)-( N T N ^ exp(-rAr(t-T))£(r)d7 



N=l 



with 



T N = 2a 1N if N + a 2N , Cn = 



2a 1N r) 2 N 



(58) 



(59) 



2a 1N i] 2 N + a 2N 

It follows from Eq. (58) that the conventional shear modulus of a rubbery polymer at small strains is 
given by 



Go 



- Ni 
PM- \ ^ PN 

15 ^ N ■ 

N=l 



This equation establishes correspondence between the material parameters at the micro-level (// and N) 
and those at the macro-level (Go and p). Formula (58) is equivalent to phenomenological relations (1) 
and (2). An important advantage of Eq. (58) compared to those equations is that the relaxation rates 
Tjy [inverse to the characteristic times tn is Eq. (2)] and the coefficients Cat [analogs of the relaxation 
strengths £jv is Eq. (2)] arc not arbitrary, but are determined by Eq. (59). 



6 Uniaxial extension of a specimen 



We consider a specimen in the form of a rectilinear bar whose points refer to Cartesian coordinates {x{\ 
with unit vectors e, . At uniaxial extension of an incompressible material, the strain tensor e is given by 



£ = £ 



e i e i ~ 2^ 6262 + 63<S3 ^ 



(60) 



where £ is the longitudinal strain. Substitution of Eq. (60) into Eq. (58) implies that the stress tensor 
reads 

a = ereiei + a {) {e 2 e 2 + e 3 e 3 ), 

where 



)-CnT n J exp(-r N (t-T))e(T)dT 
a (t) = -P(i) - ^ ]T ^ L) - OvIV /"exp(-r JV (i-T))e(r)d7 

AT— 1 L JO 



At) . _P (t) + ^£iL 

JV=1 



(61) 



It follows from the boundary condition on the lateral surface of the specimen that ao{t) = 0. Combining 
this equality with Eq. (61), we arrive at the formula for the longitudinal stress 



v(t)=Kj2jf e(t)-( N T N / exp(-r N (t-r))e(r)d7 

N=l <- ^° 



(62) 
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with _ 

k = h*. 

5 

The parameter K is proportional to Young's modulus, E 0} which is determined by the formula 

Eo = Kj2 P f- (63) 

JV=1 

To study the response of a specimen to periodic stretching, we set 

e(t) = Acxp(zwt), (64) 

where i = V^T, A is the amplitude and u> is the frequency of oscillations, substitute expression (64) into 
Eq. (62), and obtain 

Ni , - t N 

N=l I JO ) 

= A" J] 1 - Ovlw / cxp -(T N + icj)s ds U(t). (65) 
jv=i iV ^ ^0 J 



Adopting the conventional approach to the analysis of steady oscillations, we replace the upper limit of 
integration in Eq. (65) by infinity, calculate the integral, introduce the complex modulus 



and find that 



Bearing in mind that 



JVl / 

*» = *£f ( 



JVl 

Pn 1 1 c,ivi jv 



Tat + 



= E' + tE" , 

where E' is the storage modulus and E" is the loss modulus, we arrive at the formulas 

() ' () (66) 

Equations (66) are rather standard (except for the presence of the multiplier pn/N) and may be found 
in textbooks on the mechanics of polymers (see, e.g., Ferry, 1980). The novelty of our approach is that 
it provides a new scaling for the relaxation rates, Tjv, and the relaxation strengths, £aj- It follows from 
Eqs. (12), (50) and (59) that 

F N = aN 2 + bN\ Cjv = _^, (67) 

where 

a = 2 ai f = ^fj 2 , b = a 2 = (68) 

K K 

The advantage of Eqs. (14), (66) and (67) is that they contain only 5 adjustable parameters (which, in 
general, are temperature-dependent): the elastic modulus K, the average number of RHMs per strand 
in the equilibrium state N, the standard deviation of the number of RHMs S, the rate of suppression 
of RHMs a, and the rate of activation for frozen segments b (the parameter N\ is chosen in such a way 
that pNt becomes negligible compared to unity; for the quasi-Gaussian distribution (14), it suffices to set 
Ni > N + 3E). This number is essentially less than the number of adjustable constants conventionally 
employed to fit experimental data [for example, we refer to Wu and Liechti (2000), where 14 constants in 
Eq. (2), tjv, were fixed a priori and the other 14 constants, £jv, were found by the least-square algorithm]. 
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7 Comparison with experimental data 



To verify the constitutive equations (14), (66) and (67), we fit experimental data for a carbon black filled 
natural rubber vulcanizate with the filler content 50 phr (parts per hundred parts of rubber) and the 
glass transition temperature T g = 220 K. For a detailed description of specimens and the experimental 
procedure, the reader is referred to Lion (1998). 

We begin with matching experimental data in dynamic tensile tests at the temperature T = 253 K, 
the smallest amplitude of oscillations A = 0.006 and the frequency / = w/(2n) in the range from 10 _1 to 
10 2 Hz. Given a modulus K, the quantities a, b, N and S that ensure the best quality of approximation 
for the storage modulus, E', are found by the steepest-descent procedure. The parameter K is determined 
by the least-squares algorithm. Afterwards, we fix the constants a and b found by matching observations 
for E'(uj) with A = 0.006 and fit experimental data for other amplitudes of oscillations, A = 0.011, 
A = 0.028 and A = 0.056, by using three adjustable parameters: K, N and S. Experimental data and 
results of numerical analysis (with Ni = 500) are plotted in Figure 2. The same procedure of fitting 
is repeated for three other temperatures: T = 296, T = 333 and T — 373 K. Appropriate results are 
depicted in Figures 3 to 5. Figures 2 to 5 demonstrate fair agreement between observations and results 
of numerical simulation at all temperatures under consideration. 

The average number of RHMs in an equilibrated strand, N, is plotted versus the amplitude of oscil- 
lations, A, in Figures 7 to 10. Experimental data are approximated by the phcnomcnological relation 

7V = 7Vo+7ViA, (69) 

where the constants Ni are found by the least-squares technique. Figures 6 to 9 show that Eq. (69) 
correctly predicts observations in the entire rang of amplitudes, except for a region of very small strains. 
The latter may be explained by insufficient accuracy of the dynamic tests with extremely small amplitudes 
on preloaded specimens (with the compression strain of about 10 %). 

The standard deviation of the number of RHMs in a strand, S, is depicted in Figures 10 to 13 as a 
function of the amplitude of oscillations, A. Observations are approximated by the linear dependence 

S = S + Si A, (70) 

where adjustable parameters Si are determined by using the least-squares algorithm. Equation (70) 
provides an acceptable quality of matching observations in the entire region of amplitudes of oscillations, 
A, except for very low intencities of strains. 

The parameters Ni and Sj found at various temperatures are collected in Table 1 which demonstrates 
that the average number of RHMs per strand grows with the amplitude of straining (in accord with the 
intuitive picture that mechanical loading increases the molecular mobility of chains). 

The quantities a and b are plotted versus temperature T in Figure 14. Because the suppression- 
activation mechanism for RHMs is assumed to be driven by thermal fluctuations, the rates of transfor- 
mation are approximated by the Arrhenius dependences, 

a = 5cxp(-|l), 6 = 6exp(--^), 

where a and b are appropriate rates at high temperature, E a and E^ are activation energies, and R is the 
universal gas constant. It follows from these formulas that 

loga = a -^, log& = 6 -!p (71) 



where 



a = loga, ai = 6 = log&, bi = n ^ & in , log = log 10 . 



iJlnlO' ° ' i? In 10' 

Figure 14 demonstrates an acceptable agreement between experimental data and their approximation by 
Eqs. (71) with adjustable parameters en and bi, except for the region of high temperatures, where relations 
(71) overestimate the transformation rates. The latter may be explained by the following reasons: 
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1. the average radius of a tube which reflects constrains on the motion of a strand imposed by surround- 
ing macromolecules strongly increases at elevated temperatures, which implies that transformations 
of RHMs are slowed down; 

2. equilibration of an elastomer specimen at temperatures near T = 100 °C before mechanical tests 
results in a noticeable scission of chains (see, e.g., Palmas et al., 2001), which is beyond the scope 
of the present model. 

The parameter b grows with temperature faster than a (b\ = 2.32 • 10 3 versus a\ = 1.48 • 10 3 ). 
Assuming the quantity fj to depend on temperature rather weakly (as a power law, but not in the 
exponential fashion), we find from Eq. (68) that [3 substantially increases with temperature. The latter 
means that the excess free energy caused by suppression and activation of RHMs, see Eq. (30), is a 
strongly increasing function of temperature. 

Another important conclusion may be drawn from results listed in Table 1 which show that the 
parameters N\ and Si increase with temperature (the only exception is the value of N\ at the lowest 
temperature, T — 253 K, when the elastomer demonstrates quasi-brittle response). This implies that 
the influence of the amplitude of oscillations, A, on the distribution of strands with various numbers of 
RHMs becomes stronger with temperature. 

It is worth noting that the parameters N and £ corresponding to the smallest amplitude of straining, 
A = 0.006, are practically independent of temperature, see data collected in Table 2. This means that 
mechanical loading may be thought of as the main reason for changes in the distribution of RHMs. 
The effect of mechanical factors is enchanced at elevated temperatures and is weakened at relatively low 
temperatures, in agreement with the assumption about the entropic nature of transformations of RHMs. 

8 Discussion 

Our purpose now is to develop conditions under which the viscoelastic response of rubbery polymers is 
substantially weakened and to establish some correspondence between these relations and observations 
exposed in the introductory section. 

It follows from Eqs. (67) and (68) that 

Cn = (i + ^)~\ (72) 

Combining Eqs. (12), (66) and (72), replacing the parameter N by its average value, N and assuming 
Tm to be of the same order of magnitude as the frequency of oscillations (which results in a very rough 
estimate) , we conclude that the mechanical damping of an elastomer disappears (E" — > 0) , provided that 

/^ 2 (^) 2 »1- (73) 

Based on Eq. (73) one can speculate that the weak damping reported by Leblanc and Cartault (2001) is 
associated with a large average number of RHMs per strand, N, in uncured styrene-butadiene composites. 
This conclusion does not contradict the presence of filler particles in these compounds, because at the 
length scale of RHMs (which is less than the length-scale of a strand by two orders of magnitude, see 
Table 1), the influence of particles is negligible. 

Observations by Liang et al. (1999) may be explained within the framework of Eq. (73) by the large 
energy necessary for suppression of a RHM, j3, in a polymer composite with a high degree of crystallinity 
and a relatively small content of elastomer (about 10 %). 

To explain experimental results reported by Aksel and Hiibner (1996) and Clarke et al. (2000), we 
consider uniaxial relaxation test with a small strain eo, 

e(t) = (t<0), e(i) = e (t > 0), (74) 

substitute expression (74) into Eq. (62), introduce the relaxation modulus 
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and find that 



(75) 



N=l 



It follows from Eq. (75) that changes in the relaxation modulus, E(t), are negligible when either the 
parameters (,n are small compared to unity [which is described by Eq. (73)] or the characteristic times of 
relaxation tn — T^ 1 are small compared to the characteristic time of the test, r . Equations (12), (67) 
and (68) imply that the latter condition is given by 



This inequality implies that no viscoelastic response is observed an a test (even if C,n is of the order of 
unity) under the condition 



Based on Eq. (76), one can speculate that the absence of time-dependent response in experiments on 
polybutadicnc rubber at room temperature (Aksel and Hiibncr, 1996) and polyacrylate rubber at elevated 
temperatures (Clarke et al., 2000) is associated with slow rates of entropy production in these polymers 
[small ii is Eq. (48)], whereas noticeable relaxation of stresses observed in natural rubber at T — 80 °C 
(Clarke et al., 2000) may be explained by a substantial increase in the dimensionless ratio k/(fiTo) with 
temperature. 

9 Concluding remarks 

Constitutive equations have been derived for the time-dependent response of rubbery polymers at finite 
strains. A particle-reinforced elastomer is modeled as a network of long chains bridged by permanent 
junctions (chemical crosslinks, entanglements and filler particles). In agreement with the tube concept, a 
strand between two neighboring junctions is supposed to be confined to a tube which describes constrains 
on its micro-motion imposed by surrounding macromolecules. Unlike the conventional approach that 
postulates a constant cross-section of the tube, its radius is assumed to change substantially with the 
longitudinal coordinate. The wavyness of the tube's surface is associated with local heterogeneities of the 
polymer (both in space and time) which impose additional constrains on the motion of macromolecules. 
This implies that a strand may be treated as a sequence of threads (RHMs) whose thermal motion is 
unconstrained (in the sense that they go through all possible configurations during the characteristic 
time of experiment) and highly constrained segments with a negligible entropy (which correspond to the 
bottle-neck points of the tube) . Thermal fluctuations result in reconstruction of the tube surface, which 
is reflected by freezing of some RHMs and activation of frozen segments. The viscoelastic behavior of a 
rubbery polymer is associated with thermally agitated changes in the concentrations of RHMs. 

The stress-strain relations at finite strains and the kinetic equations for the numbers of RHMs in 
strands are developed by using the laws of thermodynamics. At small strains, these relations are reduced 
to the integral constitutive equation in linear viscoelasticity (1) with the relaxation kernel in the form of a 
truncated Prony series (2). Unlike the phcnomcnological approach which provides no information about 
the relaxation spectrum, a novel scaling law (67) is derived for the characteristic times of relaxation. 5 
adjustable parameters of the model are found by fitting experimental data in tensile dynamic tests for 
a carbon black filled natural rubber vulcanizate at various temperatures, amplitudes and frequencies of 
oscillations. Fair agreement is demonstrated between observations and results of numerical simulation. 
It is shown that the experimental constants are altered by temperature and amplitude of straining in a 
physically plausible way. 






(76) 
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Tabic 1: Adjustable parameters Ni and Ej at various temperatures 



T K TVp jVi S Ei 

2531) 127/77 702.07 4CL52 -16.46 

296.0 141.70 30.66 31.00 0.00 

333.0 140.00 315.69 43.58 43.58 

373.0 115.22 3139.68 25.34 600.29 



Table 2: Adjustable parameters N and E for oscillation tests with the amplitude A = 0.006 



TK N E 

253.0 122.0 33.0 

296.0 137.0 30.0 

333.0 132.0 29.0 

373.0 137.0 30.0 
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Figure 2: The storage modulus E' MPa versus the frequency of oscillations / Hz in a dynamic test with 
the strain amplitude A at T — 253 K. Circles: experimental data (Lion, 1998). Solid lines: results of 
numerical simulation. Curve 1: A = 0.006; curve 2: A = 0.011; curve 3: A = 0.028; curve 4: A = 0.056 
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Figure 3: The storage modulus E' MPa versus the frequency of oscillations / Hz in a dynamic test with 
the strain amplitude A at T — 296 K. Circles: experimental data (Lion, 1998). Solid lines: results of 
numerical simulation. Curve 1: A = 0.006; curve 2: A = 0.011; curve 3: A = 0.028; curve 4: A = 0.056 
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Figure 4: The storage modulus E' MPa versus the frequency of oscillations / Hz in a dynamic test with 
the strain amplitude A at T = 333 K. Circles: experimental data (Lion, 1998). Solid lines: results of 
numerical simulation. Curve 1: A = 0.006; curve 2: A = 0.011; curve 3: A = 0.028; curve 4: A = 0.056 
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Figure 5: The storage modulus E' MPa versus the frequency of oscillations / Hz in a dynamic test with 
the strain amplitude A at T — 373 K. Circles: experimental data (Lion, 1998). Solid lines: results of 
numerical simulation. Curve 1: A = 0.006; curve 2: A = 0.011; curve 3: A = 0.028; curve 4: A = 0.056 
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Figure 6: The average number of RHMs in a strand N versus the amplitude of oscillations A at T = 253 
K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of the experimental data by 
Eq. (69) 
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Figure 7: The average number of RHMs in a strand N versus the amplitude of oscillations A at T = 296 
K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of the experimental data by 
Eq. (69) 
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Figure 8: The average number of RHMs in a strand N versus the amplitude of oscillations A at T = 333 
K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of the experimental data by 
Eq. (69) 
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Figure 9: The average number of RHMs in a strand N versus the amplitude of oscillations A at T = 373 
K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of the experimental data by 
Eq. (69) 
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Figure 10: The standard deviation of the number of RHMs in a strand E versus the amplitude of 
oscillations A at T = 253 K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of 
the experimental data by Eq. (70) 
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Figure 11: The standard deviation of the number of RHMs in a strand S versus the amplitude of 
oscillations A at T = 296 K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of 
the experimental data by Eq. (70) 
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Figure 12: The standard deviation of the number of RHMs in a strand E versus the amplitude of 
oscillations A at T = 333 K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of 
the experimental data by Eq. (70) 



26 




Figure 13: The standard deviation of the number of RHMs in a strand S versus the amplitude of 
oscillations A at T — 373 K. Cirles: treatment of observations (Lion, 1998). Solid line: approximation of 
the experimental data by Eq. (70) 
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Figure 14: The parameters as 1 (unfilled circles) and b s 1 (filled circles) versus temperature T K. 
Cirles: treatment of observations (Lion, 1998). Solid lines: approximation of the experimental data by 
Eq. (71). Curve 1: a = 4.5289, a x = 1.4759 • 10 3 ; curve 2: b = 6.5638, h = 2.3182 • 10 3 
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